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Abstract 

Using exact diagonalization numerical methods, as well as analytical argu- 
ments, we show that for the typical electron densities in chaotic and disordered 
dots the peak spacing distribution is not bimodal, but rather Gaussian. This 
is in agreement with the experimental observations. We attribute this behav- 
ior to the tendency of an even number of electrons to gain on-site interaction 
energy by removing the spin degeneracy. Thus, the dot is predicted to show 
a non trivial electron number dependent spin polarization. Experimental test 

of this hypothesis based on the spin polarization measurements are proposed. 
PACS numbers: 73.23.Hk, 05.45. +b, 73. 20.Dx 
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The distribution of the addition spectrum of chaotic and disordered quantum dots mea- 
sured in recent experiments [|1]-|3|] seem to contradict the predictions of the orthodox constant 
interaction model The most striking feature in the addition spectrum distribution is 

the absence of any signs of a bimodal structure corresponding to the two spin states of the 
tunneling electron. In the constant interaction (CI) approximation, the ground state energy 
of a quantum dot populated by iV electrons is given by En = e 2 N 2 /2C + Ylf2i Vi + Y^i Vi 
where C is the dot's constant (or slowly varying) capacitance, r]i are the single particle 
energies, and iVj (N\) is the number of electrons for which the spin is in the up (down) 
state (N = iVf + iV[). In the absence of a magnetic field, for an even number of electrons 
jV T = JV; = N/2 while for an odd N, iV T = (N + l)/2 and iVj_ = (N — l)/2. Thus, the change 
in the chemical potential needed to add a particle, A 2 = Hn — fJ>N-i — En — 2En-i + En-2, 
is equal to e 2 /C for iV even, while A 2 = e 2 /C + V(N+i)/2 — V(N-i)/2 f° r N odd. Assuming 
that the single electron level spacings follow the random matrix theory (RMT) predictions, 
the distribution of the spacings should be, in the absence of a magnetic field, 



where s = (A 2 — e 2 /C)/2A, A = (rji — Jft-i) is the averaged single level spacing and 
(...) denotes an average over realizations or levels. As we have mentioned, none of the 
experimental distributions shows even the slightest sign of any bimodal distribution. 

When the interactions are treated beyond the CI model, different distributions emerge. 
In Ref. 0, the influence of the interactions on the distribution within the Hartree approxi- 
mation was considered. For strong interactions, this, roughly speaking, leads to a distribu- 
tion composed of two uncorrelated Wigner distributions resulting in the spin resolved (SR) 
RMT distribution (see [[||9|]) which is plotted in Fig. ^]a. This distribution begins abruptly 
(P(s < 0) = while P(s = 0) = 1/2) and is strongly asymmetric. All the experimental 
distributions are nevertheless more or less Gaussian with only weak asymmetries in the tails 



Thus, none of the above models can satisfactorily reproduce the experimental distribu- 




(1) 
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tions. We therefore use an exact diagonalization study to investigate the role of spin in 
determining the addition spectrum of a disordered tight-binding Hamiltonian with on-site 
and long range interactions. Although we pay the price of handling only small systems, we 
avoid the uncertainty of approximations and fully include the effects of electron correlations. 
In similar studies of symmetric dots, atoms and nuclei PJ7j,|ITJ[| the spatial symmetry of the 



confining potential plays an important role and leads to a strong dependence on the particle 
numbers via the shell structure (magic numbers) and a non-trivial total spin due to Hund 
rules. Because of the chaotic nature of the dots |l|-[3|] we do not expect such dependence on 
the electron number to play an important role, and therefore the study of a small number 
of electrons is still useful in understanding the properties of dots which are populated by 
an order of magnitude more electrons. Both components of the Coulomb interaction on a 
lattice, i.e., the long range part and the on-site interaction must be considered since without 
the long range component of the Coulomb interaction the classical capacitance behavior 
for the average spacing will not be reproduced, and for electrons with spin on a lattice a 
Coulomb interaction must have an on-site component. 

Three different regions of behavior emerge. The weak interaction limit, which corre- 
sponds to an on-site interaction parameter smaller than a fifth of the bandwidth, is char- 
acterized by a strong even-odd asymmetry of the addition spectrum. As discussed be- 
low this region correspond to values of r s < 0.1 — 0.3 (where r s is the ratio of the inter- 
particle Coulomb energy to the kinetic energy). An intermediate region corresponding to 
0.1 — 0.3 < r s < 1, where there is no sign of the spin in the A 2 distribution, the distribution 
tends towards Gaussian, and the total spin of an even number of electrons in the dot is likely 
to be partially polarized, while the spin of an odd number of electrons is hardly affected 
and is equal to h/2. For the strong interaction regime r s > 1 the distribution is a Gaussian 
of a width proportional to e 2 /C, and the spin polarization of the dot shows a complicated 
behavior. 

Since all experiments are in the r s ~ 1 regime these results clearly are consistent 
with no bimodal distribution being observed in the experiment. As we shall see, the unimode 



nature of the distribution can be understood by treating more clearly the ground state spin 
degeneracy. 

We calculate the ground state energies for a system of interacting electrons modeled by 
a tight-binding 2D Hamiltonian given by: 

H = e k,A,j,s a k,j,s ~ V J2 (°l,j+i,s a k,j,s + al +1JtS a kJ , a + h.c) + H int , (2) 

where e k j is the energy of a site (k,j), chosen randomly between — W/2 and W/2 with 
uniform probability V is a constant hopping matrix element, and s =f, j is the spin. The 
interaction Hamiltonian is given by: 

tiint - u 2^ a fc,i,T°^i-T a fc,ia a fc,ia + u 2^ — ,-. _ ? rjr — \A) 

k,j k,j>l,p;s,s' l' fe .J 'l,P\/ u 

where U' is the on-site interaction constant and U = e 2 /b is the long range component of 
the Coulomb interaction constant where b is a lattice constant. 

We consider a 4 x 3 dot with M — 12 sites and up to N = 9 electrons. The size of 
the many-body Hilbert space is m = thus for M — 12 and TV = 9 we end up with 

m — 1, 307, 504. One can use the fact that in the site representation the Hamiltonian matrix 
has no off-diagonal terms which couple states of different S z (where S z is the component of 
the total spin in the z direction) to diagonalize each block with a given S z separately. The 
size of each block is (jjf ) x (jlf ) which for the largest case is equal to 392,040. After the 
minimal eigenvalue for each block is obtained the global minimal eigenvalue is found. Since 
the system has a 2S* + 1 ground state degeneracy (where S is the total ground state spin), 
for all the cases in which the total ground state spin S ^ all the blocks corresponding to 
S z = —S, —S + 1 ... S have the same minimal eigenvalue. This gives an excellent check for 
the accuracy of the diagonalization procedure. Therefore, both the ground state energy and 
the value of S are found by the exact diagonalization. 

The strength of the long range component of the Coulomb interactions, U, is varied 
between — 6V. For the results presented here, an on-site coupling U' = (10/3)C7 was 
chosen, corresponding to Hubbard's calculation of the ratio of U' to U for weakly overlapping 
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hydrogen like wave-functions |TT|. As will be discussed later the exact ratio of U' to U does 
not play a crucial role in determining the main results presented here. The disorder strength 
is set to W — 8V in order to assure perfect RMT behavior for the non-interacting case. We 
also present some data pertaining to the W = 3V case, where only the lower levels (up to 
N = 6) which show RMT behavior were considered in order to gain some insight into the 
role of disorder. For each value of U, the results for each value of N are averaged over 500 
different realizations of disorder. 

The results of the tight-binding calculations at a given interaction strength U can then 
be compared to the experimental quantum dot density parameters. The ratio of the average 
inter-particle Coulomb interaction and the Fermi energy r s = Xj^fixnaB (where n is the 



electronic density and ag is the Bohr radius) corresponds to r s ~ \Jix/2(U/AV) for iV = 6, 
M = 12. For all experimental setups r s ~ 1 resulting in U ~ 4V. 

The average value (A 2 ) is presented in Fig. |l|. In the non interacting case (U = 0), 
(A 2 ) = for even electron numbers, while for odd numbers (A 2 ) = A. Even for weak 
interactions (U = 0.6V corresponding to r s ~ 0.2) there is no remnant of the even-odd 
asymmetry in (A 2 ). Moreover, up to U ~ 4, for all values of N, (A 2 ) = e 2 /C + A', where 
the capacitance e 2 /C = 0.67U was calculated through the random phase approximation 
(for details see Ref. ||12|| ) and has no adjustable parameters, and A' = 0.81A is obtained 
from a fit. Since this behavior appears already at weak interactions and the influence of the 
long range component of the Coulomb interaction is well described by the capacitance up 
to U ~ 4V (r s = 1) it is natural to concentrate on the role of the on-site interactions. In 
first order perturbation theory for the on-site interaction strength, where the single electron 
eigenf unctions are assumed to follow the random vector model (RVM), it has been shown 
|TJ that (A 2 ) = e 2 /C + 3U'/(M + 2) for an even number of electrons and (A 2 ) = e 2 /C + 
A — 2U' I (M + 2) for the odd case. Thus, the even and odd values of (A 2 ) will coincide 
at U' c = (M + 2)A/5 ~ 8V/5, which corresponds here to U — 0.5V, in good agreement 
with the numerical results. We have checked that this value of U' c holds for different ratios 
between U and U' . Thus, even for the smallest conceivable ratio (U' = U, since U' < U 



represents a locally attractive interaction) the even-odd asymmetry will disappear for r s > 
0.4, while a more reasonable estimation of the ratio will yield the values quoted in the 
introduction. Above this value of interaction it is not possible to continue to use this first 
order perturbation, but it is reasonable to assume that the excess value of (A2) above e 2 /C 
will be of order A' = 2>U'J(M + 2), i.e., A' ~ 0.6A, which is not far from the numerical 
result. Above U ~ 4V (r s = 1) short range correlations in the electronic density appear and 
the random phase approximation is no longer valid, as discussed in Refs. HTj|l2| . 

The fluctuations (5 2 A2) = ((A2) 2 ) — (A2) 2 are portrayed in the inset of Fig. [I]. As with 
the average, the even-odd asymmetry disappears at U — 0.6V and a gradual enhancement 
of the fluctuations as a function of the interaction strength is seen. The full distribution of 
A 2 obtained for all values of N (N = 4 ... 9) is shown in Fig. ^|. In Fig. [|a the distribution 
is plotted as function of A 2 — e 2 /C, which in a sense captures the distribution of the "single 
electron level" spacings. In the non-interacting case the CI+RMT distribution (Eq. |]) 
fits perfectly. Even for weak interactions, U = 0.6V (U' = 2V), the bimodal distribution 
is wiped away. Neither does the SR-RMT distribution |]3],|9| fit. The best fit for weak 
interactions is obtained from the usual RMT distribution 

Prmt(s) = (tts/2) exp(-7r S 2 /4) (4) 

where s = (A 2 — e 2 /C)/(A 2 — e 2 /C). This feature does not seem to depend much on the 
disorder, as can be seen form the W = 3V,U = 0.6V curve in Fig. Qa. As the interaction 
strength increases, the distribution becomes wider, less asymmetric, and from the fact that a 
considerable weight of the distribution is at negative values it becomes clear that no constant 
interaction model is able to describe it. For the region of r s ~ 1 (U ~ 4V) a better fit of 
the distribution to a Gaussian is obtained. The same results plotted as function of A 2 / (A 2 ) 
are shown in Fig. |2|b. A clear crossover from the RMT like distribution to a Gaussian 
distribution is seen as the interaction is increased, although its width does seem to depend 
on the disorder (see the W = 3V,U = 3.6V curve in Fig. |2|b). 

In all experiments a distribution which is Gaussian (up to deviations in the tail) is seen 
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|T|-[3|]. Since the experiments are performed in the region of r s ~ 1 this is in good agreement 
with our results. Leaving aside the width of the distribution for the moment, we need to 
understand better the unexpected disappearance of the bimodal distribution for intermediate 
coupling. We can gain insight into this by studying the ground state polarization of the dot. 
In the non-interacting case we expect the total spin of the dot to be S = for an even 
number of electrons and S = 1/2 for an odd number. For weak interactions, using first 
order perturbation in the on-site interaction strength with the RVM one finds that on the 
average the system gains 3U'/(M + 2) interaction energy by flipping one spin, while it 
loses rjN/2+i — Vn/2 kinetic energy for an even number of electrons, compared to the gain 
of 4U'/(M + 2) interaction energy and the loss of rj( N+3 y 2 — V(n-i)/2 kinetic energy for an 
odd number. The dot thus will flip a spin if the gain in interaction energy will be larger 
than the loss of kinetic energy [13| . The probability of finding two consecutive small single 



electron level spacings is much smaller than the probability of finding one small level spacing, 
therefore we expect a finite probability for an even number of electrons to be in a S = 1 state, 
while a much lower probability to find an odd number of electrons at S = 3/2 is expected. 
Indeed in Fig. |3]a it can be seen that while the average spin of the ground state for an odd 
number of electrons is (S) ~ 1/2, except for strong interactions U > AV (r s > 1), the average 
spin for the even case (S) > also for rather weak interactions. Since both the interaction 
energy and the kinetic energy scale as 1/M we expect this behavior to hold also for large 
systems. Thus, in the r s < 1 region we expect to see a significant number of S = 1 states 
for an even number of electrons and almost no S = 3/2 states for an odd number. Higher 
spin polarizations seem to be rather rare, although we encountered two realizations with 
S = 2 for N = 8 and U = 1.8V. Once interactions are stronger the even odd asymmetry 
in the ground state spin polarizability is less pronounced. This is illustrated in Fig. ^b, 
where (\Sn — Sjv— 1 1) is plotted. For any non-correlated behavior \Sn — Sjv-iI = 1/2; while 
for correlated behavior adding an electron to the dot may flip the spin of other electrons 



already in the dot and even lead to spin blockade WM. A clear indication for the appearance 
of a correlated state is seen for r, > 1. 



Recently, Stopa |L5| has suggested that spin polarization in a chaotic dot may appear 



due to scar states. According to this scenario the scar state will be populated first by an j 
electron, then several other states of higher kinetic energy will be populated in the regular 
sequence and eventually a j will populate the scar state. This scenario leads to a shift by 
a 1/2 of Sat between populating the scar state with the f electron and the j electron. This 
is a distinct spin polarization pattern from the ones previously described. The combined 
effect of scar states and regular chaotic or disordered states considered in this letter is an 
interesting question now under investigation. A very recent paper |L6] has suggested that 



the absence of the bimodal distribution is due to the deformation of the confining potential 



of the dot as electrons are added JT7[]. Again such a scenario will lead to a different spin 
polarization pattern than the previous ones, i.e., it will be equivalent to the non- interacting 
one. The measurement of the dot's spin polarization is feasible ||18||, probably via measuring 



the magnetic field dependence of the differential conductance |19|], from which the change 



in dot's spin may be observed. As shown in this letter, by studying the change in the spin 
polarization it is possible to clarify whether the underlying physics of the dot corresponds 
to the correlated regime, the intermediate regime, a deformable potential, or maybe carries 
the signature of scar states. 

There remains a question regarding the width of the distribution. While earlier exper- 
iments show a width comparable to 0.15 — 0.10e 2 /C a recent experiment || yields 



a width of 0.05e 2 /C ~ A. Our simulations show y/(5 2 A 2 ) ~ 0.2e 2 /C for W = 8V and 
(5 2 A2) ~ 0.1e 2 /C for W = 3V, consistent with previous results for spinless electron of 



(5 2 A 2 ) ~ 0.1 — 0.2e 2 /C, and with the earlier experiments (1]||. Nevertheless, for the 
numerical study at U = W (r s = 1), 0.2e 2 /C ~ A, which does not contradict Ref. [||]. It 



is clear that for r s <C 1, y(5 2 A 2 ) oc A, while for r s > 1, y (<5 2 A 2 ) oc e 2 /C. Whether the 
behavior of the fluctuations in the above model at r s ~ 1 is determined by A or e 2 /C, and 
the role of disorder, could be determined only by a careful finite size study, since the change 
in size will change the ration of A to e 2 /C and enable the determination of the relevant 
scale. 
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In conclusion, the on-site interaction is responsible of the removal of the even-odd asym- 
metry in the addition spectrum distribution. These interactions lead to a ground state 
spin polarizability of the dot. The dependence of the spin polarizability on the number of 
electrons may be used as a sensitive tool to determine the relevant physics in the dot. 

Many useful discussions on the addition spectrum of quantum dots with B. L. Altshuler, 
A. Auerbach, C. M. Marcus, A. D. Mirlin, O. Prus and U. Sivan are gratefully acknowledged. 
I would like to thank The Israel Science Foundations Centers of Excellence Program for 
financial support. 
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FIGURES 




FIG. 1. The value of (A2} as function of U for a 4 x 3 lattice and a different number of electrons 
N. All values are in units of the hopping matrix element V. Inset: the fluctuations . 
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FIG. 2. The probability distribution of A2 for different values of the interaction U, drawn on 
the scale of (a) A2 — e 2 /C and of (b) A2. The heavy lines correspond to Eq. (|I|) with e 2 /C = 0, 
SR-RMT, Eq. (0) and to a Gaussian with a variance of 0.2e 2 /C. Open symbols correspond to 
W = 8V while filled symbols to W = W. 
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FIG. 3. (a) The ground state average spin, and (b) the average of the absolute value of the 
change in the ground state spin as an electron is added to the dot, as function of the number 
electrons for different interaction strength. 
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